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Abstract. An excellent hadron to electron discrimination is a crucial aspect of calorimeter- 
based experiments in astroparticle physics. Standard discrimination techniques require full 
shower development and fine granularity but in space detectors severe limitations exist due to 
constraints on dimensions, weight and power consumption. A possible approach is to exploit 
the different neutron yield of electromagnetic and hadronic showers. NEUCAL is a light and 
compact innovative neutron detector, to be used as an auxiliary complement of electromagnetic 
calorimeters. This new approach to neutron counting relies on scintillation detectors which are 
sensitive to the moderation phase of the neutron component. The NEUCAL prototype has 
been placed after a conventional calorimeter and tested with high energy beams of pions and 
positrons. The comparison of experimental data with a detailed Geant4 simulation and the 
encouraging results obtained are presented. 



1. Introduction 

Space experiments designed to study cosmic ray fluxes at very high energies must have, among 
other features, an excellent hadron to electron discrimination. Electromagnetic calorimeters 
have a very good intrinsic hadron rejection power based on the shower topology (depth and 
lateral shape). Unfortunately, to be fully exploited, this property requires calorimeters with a 
high granularity readout and sufficient thickness to ensure the full containment of the shower. 
These features cannot be fully implemented on satellites and airborne experiments in general due 
to weight constraints and power consumption limitations. The hadron to electron discrimination 
capabilities are consequently reduced. 

A workaround solution was attempted in the successful PAMELA experiment [1]. A neutron 
detector, consisting of a moderator and a set of 3 He counters, was placed downstream of the 
calorimeter. The aim was to exploit the neutron component of the hadronic showers, much 
larger than in electromagnetic showers. This neutron counter, despite its very low efficiency, has 
been successfully used by the PAMELA collaboration for systematic effects evaluation and for 
measuring the calorimeter efficiency. 

The NEUCAL project aims to further pursue and expand the use of neutron detectors coupled 
to calorimeters by introducing the new technique of active moderation [2], that consists in 



Figure 1. CAD open view of the 
NEUCAL detector prototype. The 
scintillators are in blue, and the five 3 He 
tubes (green) are placed on top of the 
central module. 

detecting the signal of neutrons while their energy degrades within hydrogen-rich scintillators. 
The most promising time interval for the detection of the neutron moderation signal is between 
~ 10 ns to ~ 100 ns after the shower core, when the peak of the neutron flux arrives. In the 
hadronic showers neutrons within this time window are produced by nuclear excitation processes 
and have energies of the order of 1 MeV. 

2. NEUCAL prototype 

The NEUCAL prototype [3, 4], sketched in Fig. 1, consists of a three by three matrix of identical 
modules placed in a light-tight aluminum box with three shelves. Each module is made up 
of three slabs (25 cm x 8 cm x 1cm) of a fast polyvinyl-toluene scintillator (EJ-230 by Eljen 
Technology) , coupled through a common Plexiglas light guide to a fast fine-mesh photomultiplier 
(R5946 by Hamamatsu Photonics) . For the sake of flexibility, only optical grease was used in the 
optical couplings. Five 3 He tubes (12NH25/1 by Canberra) are also placed on top of the central 
module. The scintillators and the tubes were read out by fast digitizers capable of recording 
up to 10 ms of data samples. Two similar boards were used: one CAEN V1731 with 500MS/s 
capability but only 8 bit range over an input dynamics of 1 Vpp; one CAEN V1720 with a more 
limited 250MS/s capability but 12 bit range over an input dynamics of 2 Vpp. Each board has 
eight input channels. For a better timing precision one channel of each board was used to sample 
the trigger signal. Seven of the nine scintillator signals were sent to the 500MS/s board, while 
the two remaining scintillator signals and the five 3 He tubes signals were sent to the 250MS/s 
board. The readout was performed with a VME system. 

3. Testbeam setup and simulation 

In summer 2009 the NEUCAL prototype performances have been measured at CERN on the 
SPS line H4 during a test of the prototype of the the Cream-2 tungsten-scintillator calorimeter 
(CalW) [5]. Data were collected for negative pions (350 GeV), positrons (100 GeV and 150 GeV), 
and muons (150 GeV) which are not used for the present analysis. In all testbeam configurations, 
sketched in Fig. 2, NEUCAL was placed downstream CalW and all other possible devices and/or 
absorbers. In configuration a), used with positron beams, the upstream material corresponds 
to 16 Xo and 0.6 Ai; in configurations b) and c), used with pion beams, the upstream material 
corresponds to 29 Xo and lAi, and 25 Xo and 0.8 Ai, respectively. In all configurations the 
NEUCAL active volume was symmetrically positioned and orthogonal with respect to the beam 
axis. 

The NEUCAL prototype and the testbeam setup have been accurately modelled with 
Geant4 [6, 7] to validate the results against the Monte Carlo simulation. A sketch of the 
geometry of the configuration c), as implemented in Geant4, is visible in Fig. 3. The Geant4 
simulation of NEUCAL has been also cross-checked with Fluka [8, 9] with respect to the single 
neutron response, as described in [3, 4]. 
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Figure 2. Sketch of the three testbeam 
configurations used for the present analyses of 
the NEUCAL data. 




Figure 3. Testbeam geometry corresponding 
to configuration c) as implemented in Geant4. 



Testbeam data comparison with Geant4 simulations are done only at the level of the energy 
deposition in scintillators and on the number of counts in 3 He tubes. In fact, no readout chain 
response is presently modelled in the simulation. 

Two sets of simulated samples have been produced using different physics lists suitable for the 
NEUCAL use case [10]: QGSP_BERT_HP (QBERT), with Bertini model and Quark-Gluon String 
Precompound model to generate the final state for hadron inelastic scattering respectively below 
and above ~ 10 GeV; QGSP_BIN_HP (QBIC), similar to QBERT but with the Binary Cascade 
model in place of the Bertini model for the final state generation below ~ 10 GeV. Both physics 
lists include a high-precision model for low energy neutron transportation. The typical size of 
simulated samples is 20 thousand events for pions and 80 thousand events for positrons. 

4. Results 

By default data were taken in zero suppression mode, implemented in the digitizer boards, in 
order to reduce the event size and the readout time. The signals up to one millisecond after 
the initial charged particles shower were then reconstructed offline to search for delayed neutron 
interactions. The neutron signature is an isolated energy deposit in one of the scintillator 
modules or a pulse in the 3 He tubes, while a traversing charged particle will release energy on 
more than one scintillator module [4]. 

Unfortunately, the photomultipliers and the readout electronics often saturate because of the 
huge signal induced by the shower core, with different effects seen on data taken with and without 
zero suppression. In zero suppression mode the readout hardware spoilt the data in the time 
window close to the shower core as a consequence of reflections seen at the digitizer output lasting 
for a few hundred nanoseconds. In the few ten thousand events taken without zero suppression 
the board with limited ADC range (V1731) still showed slowly-recovering saturation effects. 
These issues have been understood and fixed in view of the next tests. Eventually, despite the 
difficult conditions, the samples of Table 1 have been used to produce the present results. 

In the following, all plots report points corresponding to: pion sample data (-•-) and pion 
simulation predictions for QBERT (□) and QBIC (H); positron sample data (-o-) and positron 
simulation predictions for QBERT ( ) and QBIC ( ). All plots shown in this paper are 
normalized to the number of events which develop a shower in the upstream material: almost 
all of the positron and 45% to 60% of pion events, depending on the configuration. 

The analysis of the traditional 3 He counters (see Table 1 for sample details) is an important 



Table 1. Testbeam data samples used for each analysis. ZS stands for 'zero suppression'. 



Analysis 


Time Domain 


Mode 


e + sample 


7T sample 


3 He counters 
Scint. 'early' 
Scint. 'late' 


0.1 fis - 0.1 ms 

< 1 /US 

1 /is — 1 ms 


ZS 

no ZS 
ZS 


39k, 100 GeV; conf. a) 
7k, 150 GeV; conf. a) 
39k, 100 GeV; conf. a) 


76k, 350 GeV; conf. b) 
18k, 350 GeV; conf. c) 
76k, 350 GeV; conf. b) 



validation of the simulation environment. The results are summarized in Fig. 4 where the number 
of pulses per event on all counters is plotted versus the logarithm of time, in microseconds. Both 
for data and simulation, the signal time is always referred to the arrival of the shower, taken 
as time zero. The behaviour is well reproduced by the simulation, but significant differences 
in the absolute prediction exist for pions between QBIC and QBERT, observable also in the 
scintillators analyses reported below. These discrepancies will be investigated further when new 
datasets from future tests will become available. 

The analyses of the scintillator modules data search for neutron interaction candidates (hits) 
defined as single isolated signals exceeding 0.3 GeV in energy. This definition is also implemented 
at the simulation level. All scintillator modules were previously calibrated in energy with cosmic 
muons [11]. Details of the used samples can be found in Table 1. In particular, the analysis 
of the 'early' scintillator data (< 1 fxs) is performed on samples without zero suppression; 'late' 
domain analysis, from 1 /j,s to 1 ms, relies upon data with zero suppression. 

Data from scintillators are contaminated by several background sources, not modelled in 
the simulation, for which specific rejection criteria have been deployed. Signal reflections in 
zero suppression mode are rejected by vetoing specific time delays between subsequent hits; the 
signals due to off-trigger beam particles can be reduced by asking no hits from different modules 
in coincidence; spurious effects due to saturation phenomena are reduced by rejecting signals 
with abnormally long duration. 

An example of the results of the 'late' analysis, from 1 /zs to 1 ms, is given in Fig. 5 where the 
hit energy per event is plotted as a function of the logarithm of time, in microseconds, for one 
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Figure 4. Number of pulses per event of all 
3 He counters as a function of the logarithm of 
time, in microseconds. 
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Figure 5. Hit energy per event of one 
top-side scintillator module as a function of 
the logarithm of time, in microseconds, 'late' 
domain. 




Figure 6. Hit energy per event as a function of the logarithm of time, in microseconds, for 
positrons and pions from the nine NEUCAL scintillator modules, 'early' time domain. Top row 
corresponds to the upstream module shelf with respect to the beam arrival direction. Same 
symbol conventions of Fig. 4 and Fig. 5 apply. 



top-side module; all modules show similar behaviour. This range is populated by depositions 
resulting from neutron captures on nuclei and it is not directly interesting for the detection of 
the neutron moderation. Nevertheless, the good overall agreement between data and simulation 
is an indirect indication that neutron flux estimations are under control. Capture signals could 
be used as a complementary handle for neutron detection if they can be made numerically 
significant on the single event. 

Off-trigger contribution due to beam contamination can be observed on the positron data 
for time greater than about 30 [is in Fig. 5, and also in Fig. 4. This is due to the much greater 
intensity of the positron beam used in the test. 

The 'early' time domain analysis (< 1/xs), performed on data without zero suppression, is 
represented in Fig. 6 where the hit energy per event is reported versus the logarithm of time, 
in microseconds, for all nine modules. It is evident that central modules suffer the shower core 
contamination from pions and positrons, especially at short times. Side modules are less affected 
by this issue. Nevertheless the agreement between data and simulation is again not satisfactory 
for positrons, with the exclusion of the two bottom-side modules. This is a saturation effect in 
the modules read out by the V1731 boards due to the huge signal generated by electromagnetic 
showers not fully contained by the upstream calorimeter. 

The analysis at times before 200 ns, where most of the neutron moderation signal is expected, 
is thus restricted to the bottom-side modules that are not directly flooded by the electromagnetic 
shower and that are read out by the board V1720, which has the dynamic range to withstand 
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Figure 7. Hit energy per event 
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the huge signal. The hit energy per event for the bottom-side modules as a function of time, 
now represented in linear scale, is shown in Fig. 7. The agreement with the simulation is 
satisfactory with the exception of the region below few tenth of nanoseconds where a significant 
shower contamination is still present and cannot be further reduced due to the limitations of the 
present setup. This plot demonstrates that the signal due to neutrons in pion data is significantly 
larger than in positron data in the time interval relevant for the NEUCAL application. 



5. Conclusions 

NEUCAL is a detector designed to exploit the innovative "active moderation" technique in the 
neutron counting. The final aim is the development of a light and compact device, suitable for 
space experiments, to help electromagnetic calorimeters in hadron to electron discrimination. 
The NEUCAL prototype has been tested with pion and positron showers, and data compared 
with an accurate Geant4 simulation. Results demonstrate in general good agreement between 
data and simulation and are very promising in view of the device application. Further 
confirmations are expected with upcoming tests that will profit from the experience gained 
so far and that will be organised with an improved and more performant setup. 
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